Effect of disorder on temporal fluctuations in drying induced cracking 
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We investigate by means of computer simulations the effect of structural disorder on the statistics 
of cracking for a thin layer of material under uniform and isotropic drying. For this purpose, the 
layer is discretized into a triangular lattice of springs with a slightly randomized arrangement. The 
drying process is captured by reducing the natural length of all springs by the same factor, and the 
amount of quenched disorder is controlled by varying the width £ of the distribution of the random 
breaking thresholds for the springs. Once a spring breaks, the redistribution of the load may trigger 
an avalanche of breaks, not necessarily as part of the same crack. Our computer simulations revealed 
that the system exhibits a phase transition with the amount of disorder as control parameter: at low 
disorders, the breaking process is dominated by a macroscopic crack at the beginning, and the size 
distribution of the subsequent breaking avalanches shows an exponential form. At high disorders, 
the fracturing proceeds in small-sized avalanches with an exponential distribution, generating a 
large number of micro-cracks which eventually merge and break the layer. Between both phases a 
sharp transition occurs at a critical amount of disorder £ c = 0.40 ± 0.01, where the avalanche size 
distribution becomes a power law with exponent r = 2.6 ± 0.08, in agreement with the mean-field 
value t = 5/2 of the fiber bundle model. Moreover, good quality data collapses from the finite-size 
scaling analysis show that the average value of the largest burst {A ma x) can be identified as the order 
parameter, with fi/v — 1.4 and 1/v ~ 1.0, and that the average ratio (m^/mi) of the second m^ and 
first moments m\ of the avalanche size distribution shows similar behaviour to the susceptibility of 
a continuous transition, with 7/1/ = 1., 1/v ~ 0.9. These results suggest that the disorder induced 
transition of the breakup of thin layers is analogous to a continuous phase transition. 

PACS numbers: 46.50.+a,02.50.-r,64.60.av 



I. INTRODUCTION 

Desiccation induced cracking of a thin layer of ma- 
terials is an interesting scientific problem with a broad 
spectrum of technological applications. Painted surfaces, 
thin wood layers, antirreflective and UV-protecting coat- 
ings on glasses and thin-film manufacturing processes like 
chemical-bath and sol-gel depositions are just a few of a 
large set of examples where dessication-induced cracks 
might be avoided. In contrast with fractures induced by 
collisions or tensile stresses, drying of a thin-layer ma- 
terial is a homogeneous and isotropic process occurring 
everywhere inside the layer, whose statistical behavior 
could be different from those previously studied cases. 

During the past decades experiments have revealed 
that the competition of crack formation inside the layer 
with the delamination from the substrate leads to a grad- 
ual breakup into pieces of polygonal shapes whose char- 
acteristic length scale is determined by the layer's thick- 
ness. Recently Nakahara and Matsuo addressed the pos- 
sibility of controlling the structure and time evolution of 
crack patters in pastes, i.e. dense colloidal suspensions, 
by subjecting the paste to vibration or flow before the 
onset of drying [lr[3] . Detailed experiments showed that 
those mechanical excitations induce structural rearrange- 
ments and plastic deformations getting imprinted onto 
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the patterns of the drying crack network. A large num- 
ber of experiments have been performed to understand 
the microscopic mechanism leading this memory effect. 
They conclude that the interaction (attractive or repul- 
sive) among the colloidal particles and the amount of 
disorder (shape regularity) govern the pattern formation 
[4]. A very interesting effect of the amount of disorder 
has recently been pointed out in other types of fracture 
processes, as well. Fiber bundle models with equal load 
sharing exhibit a power-law distribution of the bursts' 
sizes. 5-7 . Mixing fibers with widely different breaking 
strengths reveals a phase transition between two regimes 
with different power law exponents of the bursts' dis- 
tribution. The transition takes place at a well defined 
system composition, and it has shown to be analogous to 
a continuous phase transition |S]. 

In the present paper we focus on the effect of the 
amount of quenched disorder on the process of gradual 
breakup induced by desiccation, motivated by the works 
of Nakahara and Matsuo PUS] and Kitsunezaki (9 . The 
thin layer of paste is discretized into a regular triangular 
lattice of springs with fixed ends at the boundary. Drying 
is captured by gradually decreasing the natural length of 
all springs, inducing unbalanced forces and leading to the 
final breaking of some bonds. The quenched disorder of 
the material is represented by a random distribution of 
the bonds' breaking thresholds, and the the amount of 
disorder is given by distribution's width. There are nu- 
merous studies in the literature on the spatial arrange- 



ment of cracks [5HT2"]. Hereby, we focus on the tempo- 
ral evolution of the crack by analyzing the statistics of 
avalanches of breaking bonds. A quasi-static desiccation 
process is implemented by shrinking all the springs by the 
same rate until a first bond breaks. This induces a redis- 
tribution of loads that can be followed by an avalanche 
of breaks. Simulations revealed that the layer's breakup 
proceeds in bursts whose dynamic and statistical features 
strongly depend on the amount of disorder: At low disor- 
ders a dominating crack nucleates, creating an extended 
free surface inside the system. As a consequence, the size 
distribution of bursts is discontinuous, i.e. small sized 
avalanches have an exponential distribution while larger 
ones show a peak distribution, with a gap in between. At 
high disorders, in contrast, the breakup process consists 
of a large number of small-sized avalanches with an ex- 
ponential distribution. Varying the amount of disorder 
a transition occurs between the two phases at a criti- 
cal disorder where the size distribution becomes a power 
law. Analyzing the finite size scaling behavior of the 
largest and average burst sizes we demonstrate that the 
disorder-induced transition of our system is analogous to 
a continuous phase transitions. 



II. MODEL CONSTRUCTION 

The thin layer of material to be dried is discretized 
into a uniform two-dimensional triangular lattice with 
a slightly distorted arrangement. The nodes represent 
point masses and the bonds among them provide for 
cohesive forces. Nodes at the boundaries are set fixed, 
while internal nodes can move in the two-dimensional 
plane. System size L is defined as the number of nodes 
in the lower boundary row of the lattice. Simulations 
were carried out for the values L = 16,24,30,36,42,50. 
The bonds are simple central force springs, each one with 
a different spring constant k = El/ A, where I is the nat- 
ural length of the spring, A is a uniform cross section and 
E denotes for the Young modulus of the layer's material. 

The effect of desiccation is implemented by reducing 
the natural length of all springs by the same ratio. It is 
controlled by an increasing parameter a, given by 
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(1) 



where I is the current natural length of the spring and Iq 
denotes its initial value, computed as the initial distance 
between the two nodes connected by the spring. Springs 
are assumed to have a finite strength: When the load 
a on a spring exceeds its breaking threshold a\ h (i = 
1, . . . ,iV), the spring breaks instantaneously and it will 
never be restored (no healing). The quasi-static process 
of a very slow drying is assured by looking for the smaller 
a to break a spring. 

The quenched disorder of the material is represented by 
the randomness of the breaking thresholds a t h- In order 
to be able to vary the amount of heterogeneity in the 
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FIG. 1. (color online) Left: Example of initial randomized 
geometry. Right: Uniformly distributed breaking thresholds 
are centered at C = 1 with half- width W. The amount of 
strength disorder is controlled by varying the value of W be- 
tween and 1. 



system, we consider a uniform distribution of threshold 
values centered at C = 1 with half- width < W < 1, 
that is with probability density 
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for C - W < a th < C + W. (2) 



Hence, the strength of disorder is characterized by the 
dimensionless parameter 
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ranging from (no disorder) to 1 (highest disorder). We 
note that there is also a slight amount of structural dis- 
order in the system: Even though the spring lattice is al- 
most a triangular regular grid with spacing a, each node 
is randomly set inside a circle of radius 0.1a, centered 
at the regular place. This slight distortion is sufficient to 
prevent the formation of artificial crack patterns, but oth- 
erwise, it does not have any relevant role in the breakup 
process. 

After removing a spring, the system is relaxed by nu- 
merically integrating through the Verlet algorithm the 
equation of motion for each node, 
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where a velocity-dependent damping force (with damp- 
ing coefficient c = 10) has been added to assure for re- 
laxation. The node mass m is proportional to the area of 
the set of all points closer to this node than to the other 
ones (that is the Voronoi polygon of the dual lattice for 
this geometry). In the present paper, because the geom- 
etry is close to a regular triangular lattice, m « 0.86 for 
all internal nodes [13]. The relaxation causes the redis- 
tribution of the load dropped by the broken bond over 
the remaining ones, leading to spatial correlations in the 
system. It has to be emphasized that this force rear- 
rangement itself may result into subsequent breakings, 
as more springs surpass their breaking thresholds while 
keeping their relaxed length fixed. As a consequence, a 
single spring breaking induced by drying may trigger an 



avalanche of breakings. The avalanche size A is defined 
as the number of broken bonds during the avalanche, in- 
cluding also the first one whose breaking was provoked 
by the drying step. Thus cracks, defined as clusters of 
contiguous broken bonds, appear. Due to the long-range 
load redistribution, a single avalanche can generate sev- 
eral cracks, and it is also possible that several indepen- 
dent avalanches contribute to a given crack. 

In the framework of the model, the desiccation is 
treated as a stepwise process, where the number of 
shrinkage steps provides a measure of time t. During 
drying, the amount of shrinkage is characterized by the 
relative length reduction a = AI/Iq of spring elements, 
so that we use a as the degree of desiccation. The ini- 
tial value of a at the beginning of the loading process 
is ao = 0. The drying process ends when the humidity 
content on the sample drops to 0. We consider this to 
be equivalent to a finite value of the shrinkage parameter 
ctf, which in the present model is set to the arbitrary 
value of a/ = 0.3. A large amount of computer simula- 
tions were performed in order to understand the effect of 
the strength of disorder on the breakup process of thin 
layers of heterogeneous materials. Specifically, N = 1000 
for L = 16, L = 24 and L = 30, N = 800 for L = 36, 
N = 500 for L = 42 and N = 50 for L = 50, systems 
were generated for each value of the amount of disorder 
£ by changing the seed of the random number generator. 

In the present paper we focus on the effect of the 
amount of disorder on the statistics of avalanches. 



III. BREAKUP PROCESS 

As explained before, the desiccation-induced breakup 
of our discrete model is not a smooth process, but evolves 
in bursts (avalanches). The dessication parameter a 
grows until a single spring brakes, triggering the burst 
and leading to load redistributions and additional re- 
arrangements that eventually brake more springs. The 
burst stops when no more springs brake with that dessi- 
cation factor. Both these temporal fluctuations of the 
bursting activity and the emerging crack structure are 
strongly affected by the degree of disorder in the system. 
To give a quantitative characterization of the effect of 
disorder on the breakup process, simulations were car- 
ried out by varying the value of £ in a broad range. 

Large avalanches consist of spatially correlated bonds 
being broken; even though not necessarily as a single ex- 
panding crack of contiguous bonds. Such spatial correla- 
tion takes place at the beginning of the breaking process 
for low amounts of disorder in the distribution of break- 
ing thresholds, where large avalanches appear early on 
(as seen in Fig. [2]), due to stress accumulations. Con- 
versely, large amounts of disorder do not generate corre- 
lated crack growths during the early stages of the evo- 
lution. The initial breakings tend to be random, scat- 
tered all over the lattice, with small avalanches of size 
A ~ 1. Only after a fraction of bonds breaks and small 




FIG. 2. (color online) First large avalanche for three differ- 
ent widths of disorder, showing how the narrower breaking 
distributions favor the appearance of large fractures. Left: 
£= 0.01, Center: £= 0.1, Right:.$= 0.3. [14]. 



micro-cracks have appeared in the system, these micro- 
cracks can coalesce by means of gaping cracks, that are 
clusters of broken bonds joining two micro-cracks. As 
drying continues, both for small and large disorder, the 
newly created free surfaces hinder the appearance of large 
avalanches by relaxing the stress in the remaining intact 
parts of the spring network. The overall shape of cracks 
is similar to what has been observed for drying processes 
in composite materials, where the disorder comes from 
random distributions of the three types of bonds repre- 
senting the interaction between the ingredients of a com- 
posite (see PME2]). 




FIG. 3. (color online) The relative length reduction (a) of 
springs as a function of time for several values of the amount 
of disorder £. The system size was set to L — 30. Data values 
are binned with a bin size At — 1 and At = 16 for < t < 10 
and for t > 10, respectively. 

To follow the time evolution of the process, we evalu- 
ated the average value of the desiccation parameter (a) as 
a function of time t. It is important to emphasize that the 
curves of (a) (t) presented for several values of £ in Fig. [3] 
have two distinct regimes: at the beginning of the drying 
process the value of (a) remains nearly constant over two 
orders of magnitude in time, which implies a large num- 
ber of shrinking steps with tiny length changes 5a <C 1. 



This slow regime is then followed by a sudden acceler- 
ation of the process where springs suffer a considerable 
length reduction. The reason of the crossover from slow 
to rapid shrinkage is the formation of a dominating crack 
spanning the system in one or both directions. Due to 
the large free surface inside the lattice, a larger reduction 
in the length of all springs is required to trigger the next 
breaking event. The analysis of the spatial structure of 
damage shows that the transition point can be identified 
as the exact time when the spanning crack appears. This 
occurs for all levels of disorder £; nevertheless, the geo- 
metrical structure strongly depends on the precise value 

ofe 

It can be seen in Fig. [3] that the curves of {a} (t) shift 
to the right with increasing values of £. It means that the 
the total duration tt of the process, i.e. the total number 
of shrinkage steps required to reach some value a/, is an 
increasing function of the amount of disorder, as illus- 
trated in Fig. [4] In other words, the process slows down 



IV. AVALANCHE DYNAMICS 

Let us start by analyzing A, defined as the average size 
of avalanches decreased by the minimum avalanche size, 
A m i n = 1, as a function of a. 
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FIG. 4. (color online) Average number of steps tf needed to 
reach the shrinkage value of a/ = 0.7, for different values of 
the disorder £, for a system of size L = 30. 



with increasing disorder, i. e. a larger number of smaller 
steps are required to drive the fracture process. It can 
be observed in Fig. [4] that the monotonically increasing 
tt curve has a curvature change around £ c sa 0.5, where 
the slope of the function has a maximum. Looking at 
the rate of change, the system seems to have two phases 
depending on the amount of disorder: a rapidly evolving 
low disorder phase and a slowly proceeding high disorder 
one, separated by the inflection point. Both quantities 
we have considered until now are global characteristics 
of the lattice breaking process. In the following we ana- 
lyze the statistical features of the individual avalanches 
of spring breakings. 



FIG. 5. (color online) Arithmetic average of avalanche sizes 
A(a) as function of the shrinkage a for different values of dis- 
order f. Each curve is the average over several instances of the 
system. Inset For very high disorders (for example £= f.0), 
there is a single maximum, corresponding to a single large 
avalanche, just at the start. This first maximum disappears 
close to the transition value of disorder £ c = 0.6 and, there- 
after, a second local maximum appears. For all cases, there 
is an exponential relationship between A an a at the end of 
the process. 

It can be observed in Fig. [5] that for systems with very 
high disorder, A (a) jumps discontinuously from zero to 
an isolated maximum, which in turn is orders of mag- 
nitude larger than the average of the function over the 
whole loading process. This corresponds to a large initial 
avalanche, caused by the release of the strain accumula- 
tion in the system. The hight of this maximum decreases 
with increasing disorder and disappears for £ sa 0.7, 
where A (a) becomes a monotonically decreasing func- 
tion. This is in agreement with previous results on the 
geometrical structure of the crack patters, showing that 
for large disorders the largest crack appears by coales- 
cence, not by a single large avalanche [TUl HI] ■ Also no- 
tice that for smaller amounts of disorder (for £ < 0.6) 
a local maximum appears. This implies that the large 
initial crack does not completely undermine the appear- 
ance of large avalanches; rather, some new cracks can 
appear at this second value of shrinkage. It is interest- 
ing to note that in Fig. [5] straight lines are obtained on a 
semi-logarithmic plot, evidencing that, for large a values, 
A(a) shows an exponential decay 



A 



,-Ac 



(5) 



for all amounts of disorder. The multiplication factor A 



is practically independent of the value of £. 



a power law 
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(6) 
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FIG. 6. (color online) Avalanche size distributions iV(A) oc 
P(A) for a system of linear extension L — 30, for several 
values of disorder £. Inset avalanche size distribution at the 
critical point £ c , where a good quality power law is obtained. 

The distribution of avalanche sizes N(A), proportional 
to the probability distribution of avalanche sizes P(A), 
is presented in Fig. [6j for several values of the amount of 
disorder £. It can be observed that for small disorders 
the distribution -P(A) clearly separates into two distinct 
regimes with different functional forms: for small-sized 
avalanches (say, A < 10) the distribution P shows an 
exponential decay, while for larger ones it shows a well 
defined peak, corresponding to a single major avalanche 
at the beginning of the drying (depicted in Fig. [2]) . This 
rapid exponential decay of small-sized avalanches is due 
to the appearance of a large dominating crack of almost 
the same as the linear extension of the system, which 
creates a free surface that allows for the relaxation of 
stresses. The drying process later on cannot generate 
long breaking sequences, since nodes can easily rearrange 
to compensate for unbalanced forces. At high amount of 
disorder, in contrast, the presence of strong springs pre- 
vents the formation of a dominating crack at the onset 
of breaking. The breakup proceeds in short avalanches 
resulting in a large number of small-sized cracks which 
eventually merge and relax all unbalanced forces. Conse- 
quently, the distribution of avalanches sizes P (A) spans 
again a limited range and exhibits an exponential decay. 

These two different behaviors define the two phases 
of the system: a first isolated large crack phase where 
the fracture growth is dominated by the appearance of 
a macroscopic crack at the beginning of the process, al- 
most spanning across the sample; and a second coales- 
cence phase, where the large fracture is formed by merg- 
ing small cracks at the end of the breakup process. The 
transition between the two phases is controlled by the 
amount of disorder in the system. At the critical value 
£ c = 0.40 + ±0.01, the probability distribution becomes 



with no characteristic length scale (see the inset of Fig. 
mn. The value of the exponent r was determined nu- 
merically as r = 2.6 ± 0.08. It is interesting to note 
that this value falls very close to the mean-field expo- 
nent t = 5/2 of fiber bundles [5J Ul H], supporting the 
idea that the statistics of avalanches in drying induced 
cracking is dominated by the long range nature of the 
load redistributions. 
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FIG. 7. (color online) Arithmetic average of the largest 
avalanche, for different system sizes, as a function of the dis- 
order of the system £. The different data curves collapse in 
a single shape when the transformation is applied, implying: 
£ c = 0.4, fi/v = 1.4 and 1/v = 1. Inset Original data without 
rescaling. 

In order to understand the nature of this disorder- 
induced transition we computed the average value of the 
largest burst (A max ), which can be considered as the 
order parameter of the transition. Indeed, it can be 
observed in Fig. u\ that (A max ) is a monotonically de- 
creasing function of the amount of disorder £, with an 
inflection point at the transition (£ c rj 0.6 for the sys- 
tem size L — 48). Increasing the system size L shift 
the curves of (A max ) towards higher disorder values, but 
the functional form remains the same. Actually, Fig. [7] 
shows that all curves for different system sizes L can be 
collapsed onto a universal master curve by appropriately 
rescaling both axis on the plot. The good quality data 
collapse implies the scaling structure 
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where the value of the critical exponents were obtained 
numerically as [1/v = 1.4 and 1/v — 1.0. 

Fig. [8] compares the average of the largest (A max ) and 
second largest (A^,.) bursts. It can be observed that 
below the critical point £ < £ c the size of the largest 
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FIG. 8. (color online) Arithmetic average of the largest 
{A max } and second largest (A^^) avalanches for different 
system sizes, as a function of the disorder f. 



avalanche is a finite fraction of the total number of bonds 
in the lattice, while the second largest burst is orders 
of magnitude smaller; in contrast, in the high-disorder 
regime £ > £ c the largest and second largest bursts are 
similar and decrease together as £ increases. Based on 
this behavior, let us define the average burst size (A) as 
the average value of the ratio of the second mi and first 
moments m,\ of the avalanche size distribution 1161 , 
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where the fcth moment of the burst sizes is defined as 
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When plotted against £, this average burst size exhibits 
the remarkable feature of showing a maximum at the 
same critical disorder £ c obtained in the previous anal- 
ysis of the largest bursts (Fig. [9]) . Both results demon- 
strate that the interplay between strength disorder and 
stress inhomogeneities around cracks leads to the emer- 
gence of a critical state of the system where the bursting 
activity becomes scale-free. Note that in the evaluation 
of the average burst size (A) the largest burst is always 
omitted; thus, the maximum in Fig. [9] implies that the 
characteristic avalanche size becomes comparable to the 
system size at the critical disorder. By increasing the 
system size L the maximum of the average burst size (A) 
gets sharper, i.e. it becomes a higher and narrower peak. 
Based on the analogy to phase transitions, we tested the 
validity of the scaling form 



(A)(L,C)=i^ ff ((£-£ c )L 1 ^ ) 



(10) 



where g(x) denotes a universal scaling function. It can be 



observed in Fig. 10 that this scaling law allows for a data 



FIG. 9. (color online) Scaling plot of the average avalanche 
size (mijm\). By rescaling the two axis with appropriate 
powers of L, all curves obtained for different L values collapse 
on the top of each other for j/v — 1., 1/v = 0.9 and £ c = 0.4. 
The inset shows the original data before rescaling. 



collapse of very good quality. The value of the exponents 
were obtained as j/v = 1 and 1/v = 0.9. It has to be 
emphasized that these values of £ c and 1/v are consistent 
with the ones obtained from the scaling analysis of the 
largest burst (A max ). 



V. DISCUSSION 

We carried out a theoretical investigation of the 
desiccation-induced fracture of a thin layer of an het- 
erogeneous material. Previously, a large variety of ex- 
periments and computational modeling have shown that 
the crack patterns in such a system strongly depend on 
the amount of disorder of the material. In the present 
study we focused on the temporal evolution of the frac- 
ture process by analyzing the statistics of avalanches of 
micro-fractures varying the amount of disorder. For this 
purpose a triangular lattice of spring elements was in- 
troduced where drying was captured by gradually de- 
creasing the natural length of all springs. The quenched 
disorder of the material was represented by the random 
breaking thresholds of springs. The threshold values were 
uniformly distributed such that the amount of disorder 
could be controlled by varying the width of the distribu- 
tion. The model was investigated by carrying out exten- 
sive computer simulations for several amounts of disorder 
in a broad range of system sizes. 

As the most remarkable outcome, our study revealed 
that behind the formation of crack patterns induced by 
drying there is an interesting bursting activity of break- 
ings. These bursts are sudden avalanches of fractures, 
induced by the breaking of a single material element. 
Simulations showed that the amount of disorder has a 



very strong effect on the statistics of avalanches: At low 
disorders, few large avalanches at the beginning of the 
drying create a large free surface in the system, and there- 
after only small-sized avalanches can emerge. The large 
avalanches are comparable to the size of the system and 
they distribute like a peak around a mean size, while the 
small ones show an exponential distribution of sizes. In 
the opposite limit of high disorder the entire breakup of 
the layer proceeds in small-sized bursts with an exponen- 
tial distribution of sizes. Between the high and low disor- 
der phases a sharp transition occurs at a critical amount 
of disorder of £ c = 0.40 + ±0.01, where the avalanche 
size distribution becomes a power law. The exponent of 
this power law is r = 2.6 ± 0.08, in agreement with the 
mean-field value r = 5/2 of the fiber bundle model. 

To obtain a deeper understanding of this disorder- 
induced transition, we investigated the average size of 
the largest and second largest bursts as a function of dis- 
order and performed a finite size scaling analysis of the 
results. Actually, the average value of the largest burst 
{A max } can be identified as the order parameter, and 
a good quality finite-size scaling gives fS/v — 1.4 and 
1/v = 1.0 for the critical exponents. Similarly, the aver- 
age ratio (1712/1711) of the second 7712 and first moments 
mi of the avalanche size distribution shows a maximum 
at the same critical disorder £ c = 0.4, and the finite-size 
scaling reaches the collapse of all curves for different sys- 
tem sizes L with j/p = 1., 1/v = 0.9. These critical 
exponents characterize the transition. 

We note that in real systems the drying process might 
not be quasi-static, i.e. avalanches might be triggered 
by more than one micro-fractures, which also affects the 
statistic of avalanches. Recently, it has been shown ana- 
lytically in the framwork of fiber bundles that compared 
to the quasi-static limit, the power law exponent of the 
burst size distribution gets higher when the system is 
driven with finite load increments [6] [17]. A similar 



change of the burst size exponent can be expected when 
increasing the drying rate of the desiccation process. 

Desiccation induced cracking of pastes is a good can- 
didate for a possible experimental realization of our the- 
oretical findings. It has been shown recently in several 
experiments [5] that pastes offer several opportunities 
to control the amount of disorder in the drying material: 
By varying the composition of the system (for instance 
by mixing two different materials such as polymers and 
clay, as in Ref. [5]) a quenched structural disorder of lo- 
cal strength values can be introduced. Another possibil- 
ity is to consider a single component paste, but varying 
the size distribution of particles in the colloidal suspen- 
sion (as in [2])- Both methods might be used to prepare 
samples with several amounts of quenched disorder, and 
high-speed imaging techniques could be used to monitor 
the avalanches by analyzing the temporal fluctuations of 
the crack pattern, as in [18l|T9]. This experimental setup 
would contrast our results for the order parameter, the 
critical exponents and the phase transition itself. 
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